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A F^LY OF ^P ROX IMAT^^gOLyriONS ,^LIC TT SyROR ESTI^TES FOR THE 
NONLINEAR STATIONARY NAVIER-STOKES PROBLEM 

RALPH E. GABRIELSEN* AND STEVEN KAREL'*' 

Abstract . An algorithm for solving the nonlinear stationary Navier-Stokes 
prohlem Is developed. Explicit error estimates are given. 

1 . ^ Intrqduct i on . Since the separation prohlem of aerodynamics is at 
present intractable, it has been decided to imdertake a closely related problem 
via a mathematical technique that is potentially adaptable to the sepeuration 
problem. Specifically, the problem under consideration is the "nonlinear sta^ 
tionary Navier-Stokes problem" of fluid dynamics. The generalized Newton's 
method, as developed by Kantorovich [5, 7] is used. Its application to this 
paj^lem is of definite value for those seeking practical solutions of related 
fluid flow problems. The following qxiestions are considered: 

(i) Under what conditions does the sequence of functiohs obtained by 
Newton's method converge to the solution? 

(ii) How should the initial guess be made, as a function of v, so as to 
guarantee convergence? 

(iii) At what rate does the sequence of approximate solutions converge? 

Given S, a two-dimensional Green's domain, and fi(x,y) e CUS), 

f2(3c,y) e cMs), the nonlinear stationary Navier-Stokes problem is; 
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( 1 ) 


uUjj + + - \»Au + f i(x,y) ■ 0 , 

uv^ + Wy ♦ - V dv + f 2^*ty) “ 0 » 

X y * 

with boundary conditions 
, u(3S) « -b2(3S) , v(3S) = bi(aS) . 


An equivalent es^ression is 

(2) P(H>) = vAAip + « Alp fi - f2- ■ 0 

y A A y y a 


where 


ip^(3S) ■ bi(3S) , 41^(38) * b2( 3S) 


'I'x * V , ip^ s _u , 


and P is a mapping from C**(S) into C®(S) with norm 

3 ®ip 


UW 


N n 


na>o irnsQ 


max 


j^mg^n-m 


For clarity, the equivalence of (l) and (2) is now shown. 
LEMMA (1) (2). 

Proof. 

(1) (2) directly follows from [6, Theorem 6, p. 131]. 

(2) (l); (2) can be readily rearranged Into the form 

(-vAu + uu 4 - vu + fi) « (-vAv + uv_ + w + f2), , 

X y ^ y X y ‘ x ’ 


Let 


u ■ ip , V ■ -ip 
y * X 


z m (..vAu + UU + VU f 1 )i + (-VAv + UV + W + f2)j . 

X y X y 


Note the fact; if V = ai + bj f cMs), then 


^P^V-VF-^-t-a «b . 

y X 


with 
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Hence, q ^ ^<1. 

Therefore, (2) ■* (l). 

2 Ile Bt^t s . We seek a solution to (2) by the generalized Newton's 
method. Consider the equation 

P(t|>o) + P’ - tl'o) ■ 0 • 

If ij;i is a solution of this equation, we can write a new equation 

+ P'(rpi)(v - i|)i) = 0 . 

Assume that for each n > 0 ^ , 4" 

” n+i 


(3) 


P(ij; ) -h P'(i|» )(\p 
n n n+1 


4, ) S 0 


$ 


(See LEMMA 1 for explicit expression for P'.) 

If lia ij< exists, let be the limit. Then 
n 

n-H» 

P('}'*) + P ’ ( ) ( ^<1* - i(<*) = 0 


P(ij»*) = 0 . 


Thus ({/• is the desired solution. This is Newton's method. 
Define 


«N 



a^G(x»,y',x,y) 

3x'”3y'”-“ 


dx <3y 


where G is the Green's function of the problem 

AA$ e»0inS, iJ(3S)=0, i((^(3S) ■ 0 . 

[1]. [2], [3], [U], [8]. 

Let ■ max(max|A4>0yl . maxlAtpOxI* »axkOxl » ®«l'l' 0 y|)* 

Therefore, based on the remarkable theory developed by Kantorovich, we 
obtain the following result: 
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THEOREM 1. If the initial ie -h 




and 


as 


'^Qy 


- >>2. M H 3 < V , 
3S ^0 


|p(»o)l i 

cO 


i “♦ 0 * 3 )' 


2H 


then the Neiiiton-Kctntopovioh eequenoe fact converge to the 

unique solution of (2)» 

Proof, ly the theorem of Kantorovich [5, p. 708], it is sufficient to 
show that 

|P(r|>o)fl i 5 i 

2 |P’(\|»o)"^l ip'i 

this is shovm by the following lemmas. 

LEMMA 1. P'(ij>), the Freohet derivative exiete at all pointe in the 

domain, and 

P'(t|'o)e = + '^Oy^^x “ '^°x^*y • 

Proof: 

P(ij;) “ ’ 

Urn |?(r») - p(») - M , 0 ^ 

|a|^ 

for some linear operator L, then define P’(i|»o) * 

p(\j> 0 +a) - P{<|)o) = vAda + 'I'Oy^ftjj + ^y^'^Ox "*■ ~ ®'x^'*'°y " ’^'x^y ” “'x^y * 

Let 

L(\|»o)(a) « vAAa + 'I'OyAa^ ®'y^'l'0x " “x^'^Oy " '*'®x^y * 

Then L(ij<o) is a linear operator. Therefore, L(\|)) * P'(ip) if 


lim 

M-K) 


|p(tPo+ft) - P(’l'o) - 


|a| 


la da - a Aa J 

u„ ' y , y , . - yl' . . 0 . 

1 . 1-0 


a Aa^ - a„Aa e C° . 
y X X y 
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Therefore, 


I a Aa - a Aa I ■ max a Aa - a Aa 
yx xy" glyx xy 


a e 


C“ . 


Therefore , 


4 n 

EE 

nao n^O 


max 


S“a 


3*“3y”-‘ 


V*yl ^ I'VKI'xxxl * l“■xyy'’ * * l”-yyyl> ' 


Therefore, 

max I a„Aa -a„Aa„ I < max j a„ I f max I a | -Haax | a, 


y X X y 


(” 


s y'Vs ' 


l^+m^|ajjl^i 


max a -H&ax a 
S ' xxy' g ' - 


yyy^) 


So: 


0 < lim 
!lal!-K) 


< lim 
(afl-K) 


|p(ij> 0 +a) - P(t(»o) - La 


maxla 1 

(maxla | 

\ 

max|a | 

* maxja | 

(max|a _ | 

+ maxla 1 j 

s y 

\s 

S 

s * 

Vs ^ 

s ^ 1 


< lim Hai(lall * Mall) |ag(jal * jail) 
|a|>0 


Therefore , 

llPCifo 4- a) - P(i;»o) - La|| 
lim ' a 0 . Q.E.D. 

Ia«^ 

LEMMA 2. P”(i|^) exiete at all points c‘*(s) and 


Proof. By definition, P’*('i»o) exists if there is a bilinear operator B 4* 


lim 


|P’(K »0 + ♦) - P'('i»o) - B*l| 

Ri 


0 . 
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If ao, P"(t|»o) is defined to be B, 

[P'(t|» 0 -^^) - P’(«l»o)]« - ♦yAgx + «yd* 3 j - gjjMy - (►/gy . 

Now, 

[P' ('l>0+^) - P' (4»o)3(gi+g2) ■ [P*(t|>0+#) - P'(’Po)3gl + [P'Ci^oH) - P'('l'0^3g2 » 

and 

[P»(ij»O+<>l+02) - P'(tl»o)3g = [P’(’J’ 0 +'l>l) - P’(’Po)3g + [P'('l» 0 +'l* 2 ) - P'('l'o)]g • 

Therefore, [P*(ij>o+^) - P’(»|^o)3g is a bilinear operator of (|> and g. Let 
B4> = P'(iJ>o+<>) - P'(<Po3* Th®n i"*- follows that 

LEMMA 3. P'C'i'o)’"^ exieta. 

Proof. Given P(ij/o) + P* (i^o) “ 0, let ^ then 

P'('j»o)i’ “ -P('l>o) • 

Equivalently, vAA^» + Ai/fo^j^y + 'i'Oy^^'x " ’^'Ox^^y " 
can be abbreviated as vAAiJi ■ f(iji) + F. Let G 
ij^ (see Theorem l), then 

Ss-fGf + ifGF 

V * V » 

Define the linear operators 

Aij; » ^ J Gf(5») . 

B[-P(il;o)] =* “ • 

Then 

- ■” A^i< a B[-P(»|)o)] • 

Under proper conditions, as shown later in the proof, A^ exists. Then 


Ijj 0 -P(\|;o), This equation 

y ^ 

be the Green's function for 
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Therefore, 


Sji • (l - ^a)" B[-P(i|)o)3 . 

P’(«o)(l - ^A)'^B[-P(t|»o)3 - -P(i|>o) ; 

a)"^b - I , 

and (l - ^ a) ^B Is a right inverse. Also, (l - ” a) BP*(^q)$ ■ 

i A^*^BP' (t^o3 ■ ^ ^ inverse. Therefore 

(l - ^ a)’^B - [P’{t|^o)3“^ . 

We now show the conditions under which ^ A^ exists. 

All) « J G(i|)^hi|)0y - 


I Alj » sup 

m<i 


so 


|a| = sup f max f G(i|) ht|)(» -\|) ht|)o +)|)o„AiP -^I'Cvdil) )dS + max ( ) , 

C JiHlilLx'.y’eS'^S * xy-y x 

|1A| < sup max f [ IgI + |g , i + 1 g , | -t- . . . 3 • [ 1 1 I + • 

c ,y'ES^ ^ 

MaI < max ( max M f [ |g| '• |g , | + |g , 1 + . . .3 • [!’!' I + N 
C W<l\x',y'eS X y x y 


Therefore , 


IIAII 


/ [|g1 + Ig^.I 4 |Gy,| ♦ . . .3dx dyj . 


lAt ^ M, H, . 
C ** '*'0 


Therefore, (l - “ a)*^ exists if < v. (l - ^ a)"^ 


nsQ 


4 . . . 


. .3 


4 . . . 


exists. 
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LEM4A k, )P"(t|»o)| < 1. 

Proof, 

I!p"(*o)| - "WX * K^K - M®v - 

1^ *.^y* y* *y *y 

mh 

S ®ax(|^y| + l^jcxx^ **' l^xyyl ^ * l^xxy^ ^ ^ ^ * 


This cQo^letes the proof of Theorem 1. 

THEOMM 2. Under the l^potheeie of Theorem 1, the error eetimate for the 
mth approximate solution i,e expreeeed by 


- y 




ipi^n 


2ffl 


Proof, 5y the Kantorovich theory, this result follows from the l^vpothesis 
of Theorem 1. 


COROLLARY 1. 


Ae a function of v, for fixed m. 



COROLLARY 2. If the hypotheaie of Theorem 1 ie eatiafied, then for a given 
t > 0, there exiata a denumerably infinite number of linear equationa and aolu- 
tiona aa apeoified by (3) euoh that the entire family of ili^'e are within 
the t-neighborhood of the exact solution ^ of (2), i,e,. 


!!♦ - *„# 


c** 


< E 


Proof, Follows directly from Theorem 2, 
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